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We study the antiferromagnetic spin- 1/2 Heisenberg model on a two-dimensional bipartite quasiperi- 
odic structure, the octagonal tiling - the aperiodic equivalent of the square lattice for periodic sys- 
tems. An approximate block spin renormalization scheme is described for this problem. The ground 
state energy and local staggered magnetizations for this system are calculated, and compared with 
the results of a recent Quantum Monte Carlo calculation for the tiling. It is conjectured that the 
ground state energy is exactly equal to that of the quantum antiferromagnet on the square lattice. 



PACS numbers: 75.10.Jm, 71. 23. Ft, 71.27.+a 

In this paper a renormalization group transformation 
is used to study the ground state of Heisenberg spins 
with antiferromagnetic couplings on a two-dimensional 
quasiperiodic tiling. This system poses a novel theoreti- 
cal problem, namely, the nature of quantum fluctuations 
in a structure possessing a number of exact symmetries 
but no translational invariance. While periodic systems 
and disordered variants thereof have received much at- 
tention, little is known about aperiodic quantum mod- 
els in two or more dimensions. In particular, the real 
space magnetic ordering of local moments in systems with 
quasiperiodic long range order remains to be elucidated, 
and should present novel and complex features, different 
from properties of crystalline or disordered systems. The 
archetypal nonfrustrated two-dimensional antiferromag- 
netic system is that of spins on the square lattice, an old 
and until recently controversial problem, while the prob- 
lem we consider now , with its fundamentally different 
symmetry properties, aims to understand a new class of 
unfrustrated systems. 

Experimental work providing motivation for the study 
of such systems comes from neutron scattering studies of 
the magnetic phase in a Zn-Mg-Ho quasicrystal [1]. The 
magnetic diffuse scattering of the low temperature phase 
shows an icosahedral symmetry, reflecting the underly- 
ing quasiperiodicity of this compound. The nature of 
the ground state in such a quasicrystalline medium was 
recently discussed in [2] where Quantum Monte Carlo 
(QMC) calculations were carried out for an antiferro- 
magnetic Heisenberg model on one of the simplest two- 
dimensional quasiperiodic tilings available, the octagonal 
tiling. This tiling has been frequently used for numerical 
investigations of the effects of quasiperiodic modulations 
in two dimensions. More detailed, analytic and numeri- 
cal results are available for one dimensional quasiperiodic 
models, where quantum spins have been considered using 
real space renormalization transformation [3] , and using 
density matrix renormalization or by using mappings to 
fermionic models (see [4] and references therein). How- 
ever, the techniques used are particular to one dimen- 
sion and not readily generalisable to the two-dimensional 
structure considered here. 



The model considered in [2] has a Hamiltonian H = 
JYlu,j) Si.Sj where the spins are located on vertices of 
the octagonal tiling, and J is coupling along each edge. 
Spin-spin correlations in the ground state were com- 
puted, and the staggered local moment at a given site 
was found to depend on the number of nearest neighbors 
z. We recall that the octagonal tiling has six z-values 
ranging from 3 to 8 (see ref. [2] for a picture). Fig.l shows 
the six types of nearest neighbor configurations that oc- 
cur, along with the nomenclature used in this paper (see 
below) . 




FIG. 1. The six nearest neighbor environments on the oc- 
tagonal tiling 

Within each of the six families the local order param- 
eters have a fine structure reflecting the differences in 
the next-nearest neighbor shells, and there are further 
splittings due to even longer range spin-spin interactions. 
The observed z dependence of the local order param- 
eters was qualitatively explained in [2] by an isolated 
cluster or Heisenberg star approximation. Doing bet- 
ter requires taking into account successive shells of next- 
nearest neighbors. The best way to do this is by a renor- 
malization group calculation, using the self-similarity of 
the quasiperiodic structure. This is the aim of the present 
calculation. 

The renormalization approach is a natural one for 
self-similar quasiperiodic tilings invariant under a scale 
transformation or inflation such as the one described by 
Gardner [5] for the Penrose tiling. For the octagonal 
( Ammann-Beenkkcr) tiling [6] , one can start with tiles of 
some given edge length, and reconnect a certain subset 
of vertices (inflation). The redrawn tiling is then sim- 
ilar to the original one, except for an overall scale fac- 
tor, equal to the golden mean r = (\/5 + l)/2 for the 
Fibonacci chain, 2d Penrose tiling and its 3d general- 
ization, or the silver mean A = (1 + \/2) in the case of 
the octagonal tiling. This structural property of tilings 
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has been often exploited in order to establish recurrence 
relations for parameters occurring in discrete spin mod- 
els, electron hopping models, etc, as mentioned before for 
the one-dimensional case, and for some two-dimensional 
models [7,8], where analytical methods remain hard to 
implement. 

Our renormalization group is a generalization of the 
calculation of Sierra and Martin-Delgado for the square 
lattice [9] , where the authors considered blocks composed 
of five-spin star-shaped clusters. On the quasiperiodic 
tiling, the choice of block spins is suggested by site behav- 
ior under an inflation operation. Inflation results in the 
disappearance of low z sites. After inflation the sites that 
remain are those of high coordination number, A,B,C,and 
D\ sites, having z = 8, 7, 6, 5 respectively. We refer to 
them collectively as a sites. Sites that disappear have 
z = 5,4 and 3 respectively (the -D 2 ,E and F sites). Note 
that there are two types of five-fold site, that behave dif- 
ferently under inflation [10]. The relative number of sites 
of each kind, /j is preserved under inflation, whereas the 
density of sites is reduced by A~ 2 . After inflation, sites 
have new coordination numbers z 1 as indicated in the 
following list of transformations : 

A — > A or B or C or Di ; 

B^L> 2 ; C -> E ;L>!^F 

A natural choice for block spins is to consider star- 
clusters centered on sites of the a class. After inflation, 
the old z-blocks will become the vertices of the inflated 
tiling, and new blocks defined at the high-z sites, and so 
on. The block spins and the couplings will renormalize 
to site-dependent values, grouped according to the lo- 
cal environments. Fig. 2a shows a central Z?i site, which 
transforms after inflation to a z 1 = 3 site. The sites re- 
maining after inflation are shown with large dots, and 
the dashed grey lines represent effective interactions be- 
tween these sites. Intra-block couplings are shown by 
thick lines, while inter-block couplings are shown by thin 
dotted lines. 




FIG. 2. Block spin centers (filled circles) showing the 
central and all peripheral blocks for three cases: (a) a 
z — 5, z = 3 site (b) a z = 6, z' = 4 site (c) a z = z' = 8 site 



For an isolated block spin with z spins surrounding a 
central spin and antiferromagnetic interactions, the clus- 
ter has a spin of S' = (z — 1)5 in the ground state. The 
energy of the isolated block can be exactly found, while 
for the inter-block couplings, we will follow the approach 
used in [9] for the square lattice (where all blocks carry 



the same value of z = 4) , where one finds 5' = 35. The 
spin renormalization factors are taken to be equal to the 
classical value ^ — 1/(2—1) for simplicity. The new 
block spins 5' are situated on the black circles repre- 
senting the sites of the inflated lattice, while all of the 
nearest neighbors are decimated in the renormalization 
group (RG) transformation. 

However, on the octagonal tiling, blocks are not all iso- 
lated or disjoint as Figs. 2b and c show. Blocks as defined 
above can share pairs of sites, with a finite frequency of 
occurrence. This makes the tiling, a true two-dimensional 
structure, harder to solve than a fractal Sierpinski-type 
structure, which would have less connectivity. Fig. 3 
shows the tiling with grey dashes connecting such pairs 
of shared sites. To disconnect the clusters along the grey 
lines, the two spins are assigned to one or other of the 
overlapping blocks. This is done by annulling one of the 
couplings for each of the spins. Thus we have a diluted 
version of the original tiling, with certain couplings an- 
nulled. The fraction of annulled couplings is finite, and 
can be calculated exactly to be V2/X 3 w 0.10 or ten 
percent. 




FIG. 3. Tiling showing block centers (fat dots). The grey 
lines connect pairs of sites that are shared between two blocks. 

For example, Fig. 2c shows eight overlapping D\- 
clusters surrounding an A site. They are decoupled by 
annulling one of the couplings on either side of each D 
site. The result is that the coordination numbers of all 
Di sites after dilution goes from z = 5 to z — 3. Simi- 
larly, C sites have their z reduced from 6 to 5. 

For the seven local environments on the octagonal 
tiling, one has specific connectivity rules: an A site is 
always coupled to eight F sites, a B site always coupled 
to five F sites and 2 E sites, etc. The new block spin 
variables take on environment dependent values, and af- 
ter one inflation, one finds that = (S^\...,Sp^) — 
CS(°), with 



/ 


-1 

















8 \ 




-1 

















8 




-1 

















8 




-1 

















8 







-1 











2 


5 










-1 








3 


2 


V 











-1 


1 


2 


0/ 



(1) 



where S^ ' = (s , s , s ) and s = i. One can follow 
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the renormalisations as the tiling undergoes successive 
inflations, with = CS^™ -1 ), and it is easy to show 
that for large n w SS^" 1 ), and spins tend to relative 
asymptotic values given by(l,l,l,l,l,|,|). In a star 
cluster where all spins have different lengths, with z spins 
of Si = riiSo are coupled by the same J to a central spin 
So = «o so the ground state energy is (for n i > n o) 



e(J,z,{n}) = ~n a J(^2', 



2)/4 



(2) 



If n = (r»A, TiB,nc,nm) are the number of blocks in a 
given region of each given type, the number of blocks of 
each type after one deflation is iVn where 
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(3) 



whose largest eigenvalue is equal to 7 so that the total 
number of blocks increases (decreases) with the number 
m of deflations (inflations) as 7 m for large m. 

The effective interaction between block spins is deter- 
mined by inspecting how links transform under infla- 
tion. A minimal model can be defined by considering 
just five types of links. The set of couplings retained in 
the model is represented in an array j = (j a F,jaE,jD 1 D 2 i 
jD 2 F Jef)- Here, j a p is used to denote the link between 
(A,F), (B,F),(C,F) and (£>i,F) pairs. Similarly, j aE de- 
notes the link connecting (B,E),(C,E) and (D\,E) pairs. 
After inflation, the new couplings between sites are writ- 
ten in terms of the five old couplings, giving rise to a mul- 
tiplicative renormalization scheme [11]. After one step of 
inflation the new couplings (the grey lines in Figs. 2) are 
found to be jW = M^jW where 
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with the initial condition (taking the zero order coupling 
J = 1) = (1; f j f 1 1 • 1)- After one inflation the new 
Hamiltonian is written using averaged values of the renor- 
malization factors ^ and couplings since the effective 
spins and couplings arc no longer uniform. Block ener- 
gies are obtained using Eq.2 with a site-averaged value 
of the coupling, calculated appropriately for each of the 
seven families of sites (neglecting minor differences of lo- 
cal environment in some cases). A sites, for example, 
have eight A-F links to their neighbors, so their average 
local coupling is — j^p. The matrix evolves 
under successive inflations to a fixed point whose maxi- 
mum eigenvalue 75 w 0.15. Thus for large n, couplings 



decay as j( n ) = 75j^™ -1 ', while the corresponding eigen- 
vector determines the fixed point relative couplings. 

With these definitions, we now turn to the results ob- 
tained. The first quantity of interest is the ground state 
energy per site, eo- The QMC data in [2] yield a value 
of e ~ —0.66, while that of the square lattice was deter- 
mined numerically [13] to be about —0.67. A plausible 
conjecture is that the octagonal tiling, with its two sub- 
lattice structure and its average coordination number of 
4 has the same GS energy as the square lattice, however 
this remains to be proven. In the RG scheme, the ground 
state energy can be written as an infinite sum 
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using the fact that the frequency of blocks of type i is 
initially fi [12], and is diminished by 1/A 2 at each step of 
RG. Here e(") is shorthand notation for the energy of an 
nth stage block having a spin Sq"^ at the center, an aver- 
aged coupling value , and surrounding spins of value 
Si as given by Eq.2. The series for the energy gives 
e ~ — 0.5f (compared with the result of about -0.54 [9] 
for the square lattice). One reason for the discrepancy 
between renormalization and Quantum Monte Carlo re- 
sults for the quasiperiodic model is the appreciable bond 
dilution occurring at C and D sites, which leads to hav- 
ing fewer energy terms in the Hamiltonian. On the other 
hand, the loss of bonds is partly offset by the fact that 
the dilution tends also to suppress frustration and raise 
the local order parameter. A crude way to put back the 
"missing bond-energies" is to add in half of the miss- 
ing link energies at each of the C and D sites. This is 
easily done by adjusting the z values (zq goes up from 
5 to 5.5 while Zdi is increased from 3 to 4). One then 
finds an adjusted ground state energy of about —0.59. 
This correction technique will be applied to the calcula- 
tion of local order parameters discussed below, with good 
results. 

The QMC data in [2] give values of local order pa- 
rameters defined in terms of the local energies Ei = 
J "l2(Si-Si+s) i where the sum is over all nearest neigh- 
bors of a given site i and the spin correlations are evalu- 
ated in the ground state. The relation taken for the local 
order parameters is m s j = Ei/ z [14]. We are therefore 
interested in the cluster energies, as a function of z. 
fn the zeroth approximation the cluster energies are just 
the HS e|°^. The RG allows us to calculate the cluster 
energy 

E in) 

on increasingly bigger length scales, where z 
is the coordination number of the central spin at the end 
of n steps. For example, taking n = 1, consider an A-site 
of the inflated tiling with eight F sites around it. The an- 
cestor of the central A-site is an A site with an associated 
block energy of e^. The ancestors of the neighbors are 
D\ sites, and each one contributes half its block energy 
to the cluster energy. The total energy E^ is therefore 
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the sum of the Heisenberg star energy for A sites having 
a first-order coupling plus a zeroth order A block 
energy term, plus half the zeroth order block energies of 
its neighbors as follows: 



(6) 



and similar expressions are written for the other six types 
of site. At each stage of RG, the cluster energies are 
used to find the corresponding values of the local order 
parameters. For n = 2, cluster energies for the twice- 
inflated tiling can be written out in terms of the energies 
e\ k ^ (k — 0,1,2). The number of terms contributing to 
the cluster energy is governed by the largest eigenvalue 
of N, so that EM /7 n tends to a limit as n — > oo. In 
that limit, the quantities m s = y/ r^EW /7™ therefore 
have asymptotic values which are compared to the avail- 
able numerical data. In fig. 4a wc have compared the 
m a obtained after zero (dashed line), one and two and 
three RG steps (open circles, squares and filled circles 
respectively). After two steps, the values of m s con- 
verge quickly. In Fig. 4b are shown the third (circles) 
and fourth order (squares) results which overlap on the 
scale of the figure. The limiting values of m s are clearly 
below the QMC data, shown as grey circles, and this is 
expected due to the bond dilution. If the bond reduction 
is compensated by putting back half the bonds as we did 
earlier to estimate eo, we get estimates for m s values on 
the original octagonal tiling. The grey squares of Fig. 4b 
were obtained by correcting the n = 4 data in this way, 
in fairly good agreement with the QMC data. 
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FIG. 4. m s values plotted versus z for increasing orders of 
RG. (a) zero (dashed line), 1st (open circles), 2nd (rectangles) 
and 3rd (filled circles) order RG. (b) QMC data (grey circles), 
3rd (circles) and 4th (rectangles) order RG, and corrected 4th 
order data (grey rectangles). 



In conclusion, we have presented an RG scheme for 
a two-dimensional quasiperiodic tiling that can be com- 
pletely solved, under certain approximations. The prob- 
lem is of importance as being the simplest 2d aperiodic 
quantum antiferromagnetic spin model possible, like the 
square lattice antiferromagnet for periodic systems. The 
results of the approximate RG as regards the local or- 
der parameters are close to those calculated for the full 
undiluted model, and we believe the model takes into ac- 
count important aspects of the quasiperiodic geometry 
of the tiling. The method is less good at obtaining the 
ground state energy. We note that the ground states of 
the octagonal tiling and the square lattice appear to share 
the same value of the ground state energy - a conjecture 
awaiting proof. One notes finally that the RG calcula- 
tion after appropriate modifications can be extended to 
electronic, vibrational and other discrete problems (work 
in progress). 
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